Charge frustration and quantum criticality for strongly correlated fermions 
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We study a model of strongly correlated electrons on the square lattice which exhibits charge 
frustration and quantum critical behavior. The potential is tuned to make the interactions super- 
symmetric. We establish a rigorous mathematical result which relates quantum ground states to 
certain tiling configurations on the square lattice. For periodic boundary conditions this relation 
implies that the number of ground states grows exponentially with the linear dimensions of the 
system. We present substantial analytic and numerical evidence that for open boundary conditions 
the system has gapless edge modes. 



Quantum criticality is a key notion in the analysis 
of non-Fermi liquid behavior of strongly correlated elec- 
trons. In general, quantum criticality is expected when 
parameters are tuned such that competing orders are per- 
fectly balanced In two or three spatial dimensions, 
the identification of quantum critical behavior is quite 
difficult, due to a lack of sufficiently powerful analytic 
methods and to the limited reach of numerics. 

In this Letter, we identify quantum critical behavior in 
a specific model for spinless itinerant fermions on a two- 
dimensional square lattice. The Hamiltonian combines 
standard hopping terms with strong repulsive interac- 
tions. At filling / = 1/2 (one fermion per two sites) this 
model develops a Mott insulating phase with a fermion 
on every other lattice site. In 0], it was shown that a re- 
lated model exhibits a stripe phase upon modest doping. 
Here we study the model in the 'strongly doped Mott 
insulator phase' at intermediate density. 

In the model at hand, interactions have been tuned to 
obtain a property called supersymmetry, which we em- 
ployed in our analysis. In previous work (for a review see 
[3]), we demonstrated that this model develops a form 
of quantum charge frustration that leads to an extensive 
ground state entropy [1, Q . Here we show that this same 
model on the square lattice with an edge possesses quan- 
tum critical edge excitations. We find that this quantum 
criticality arises due to the competition of two forms of 
charge order at densities f ~ 1/5 and / = 1/4. 

Our degrees of freedom are spinless fermions living on 
the square lattice. A fermion at site i is created by the 
operator c] with {ci, cj} = Sij. The fermions have a hard 
core, meaning that they are not only forbidden to be on 
the same site as required by Fermi statistics, but arc also 
forbidden to be on adjacent sites. Their creation operator 
is dl = clv^iy, where 

r<^>^ n (1-4^^) 

j next to i 

is zero if any site next to i is occupied. 

Our model, first introduced in [6|, is easiest to define 
in terms of the "supersymmetry" operator Q = X^i*^!- 



The Hamiltonian is 

H = {Q,Q^= J2 4d,+Y.'^<^>- 

<ij> i 

The model is supersymmetric because = (Q^)^ = 0, 
which then implies that [Q, H] = [Q\ H] = 0. The latter 
term in the Hamiltonian combines a chemical potential 
and a repulsive next-nearest-neighbor potential. On the 
square lattice it reads 

i i 

where V^<i> -I- 1 is the number of particles adjacent to i, 
unless there are none, in which case V<j;> =0. is the 
number of sites. The operator F = d\di counts the 
number of fermions and commutes with the Hamiltonian. 
The chemical potential is fixed at = 4 and there is no 
a priori constraint on the particle number. 

The problem of counting the ground states in this su- 
persymmetric lattice model turns out to be related to 
some simple-to-describe (but often difficult to solve) ge- 
ometrical counting problems (see 0, A heuristic way 
of understanding this is from the "3-rule" : to minimize 
the energy, fermions prefer to be mostly 3 sites apart 
(with details depending on the lattice) . Using techniques 
from cohomology, it was proved for several lattices that 
the configurations which satisfy the 3-rule are in one-to- 
one correspondence with ground states Q • 

The correspondence between these configurations and 
ground states holds for a wide variety of boundary con- 
ditions. This indicates that the ground state wavefunc- 
tions exhibit a form of charge ordering. This ordering 
is akin to the Neel ordering in an antiferromagnet: the 
ordered state is not an exact eigenstate, but one can find 
a non-vanishing antifcrromagnctic order parameter for 
the ground state. In this sense, the ordered state dom- 
inates the ground state. For our model, the structure 
is even richer, because there are typically an exponen- 
tially growing number of ways of satisfying the 3-rule. 
Such behavior is what we mean by charge frustration. On 
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generic lattices, the number of ground states grows ex- 
ponentiallywith the fuU (two-dimensional) volume of the 
system 0,0], see also 0- However, for the square lattice 
described here, we will show that for periodic boundary 
conditions the number of ground states grows exponen- 
tially with the linear dimensions of the system. 

An important object in supersymmetric theories is the 
Witten index W, which is defined as VF = N), — Nj, 
where Ni, is the number of bosonic ground states (those 
with an even number of fcrmions), and Nf is the number 
of fermionic ground states (those with an odd number of 
fermions) Q. By definition, is a lower bound on the 
number of ground states. 
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FIG. 1: The diamonds on the left and squares in the mid- 
dle. On the right we show the square lattice rotated by 45°. 
Sublattice Si is indicated by the circles and sublattice S2 is 
indicated by the drawn lines. 



We show in this paper that for the square lattice, the 
ground states are in one-to-one correspondence with con- 
figurations that come from tiling the plane using the four 
rhombuses illustrated in fig. [l] which for obvious reasons 
we call diamonds and squares. Here the particles sit at 
the corners of the rhombuses, along the edges they are 3 
sites apart. A key result was proved rigorously by Jons- 
son 0]. He showed that for the square lattice with peri- 
odic boundary conditions in both directions, the Witten 
index is simply related to tilings with these rhombuses. 
Precisely, let tb [tf) be the number of ways of tiling the 
torus with these four rhombuses, so that there are an 
even (odd) number of fcrmions [or equivalently, an even 
(odd) number of rhombuses] . The most general version of 
this theorem allows for any type of torus on the square 
lattice with periodicities u = (mi,U2) and v = (vi,V2)- 
The expression for the Witten index then reads [9[ 

Wu,v ^Nb-Nj = u-tf- i-if-ed_ed+, 

where d± = gcd(ui ± U2, fi ± V2) and 

2 if c? = 3fc, with k integer 
— 1 otherwise. 

A natural extension of Jonsson's theorem is to relate 
the total number of ground states to rhombus tilings: 



Nb + Nf=tb + tf + A, 



where |A| = \9d_ed 



When 



u 



(1) 

(m, — m) and 



vi + V2 = 3p we prove this relation explicitly with 
A = -(-l)(^"+i'P6'd_6irf^ (see below). 



In our supersymmetric models, the problem of com- 
puting the total number of ground states reduces to find- 
ing the cohomology of the supercharge Q @ . More pre- 
cisely, the dimension of the non-trivial cohomology of Q 
corresponds to the total number of linearly independent 
ground states. The cohomology is the vector space of 
states which are annihilated by Q but which are not Q 
of something else. 

To compute the cohomology is in general very difficult. 
A useful theorem is the 'tic-tac-toe' lemma of [l^. This 
says that under certain conditions, the cohomology Hq 
for Q ^ Qi + Q2 is the same as the cohomology of Qi 
acting on the cohomology of Q2- In an equation, Hq = 
Hq^{Hq^) = H12, where Qi and Q2 act on different 
sublattices Si and 82- 

The crucial step is to choose the right sublattices: for 
the square lattice we take a set of disconnected sites for 

51 and a set of (periodic) chains for ^2 (see fig. [J). For 
free boundary conditions in either one or both of the diag- 
onal directions [(m, — m) and (n, n)] the full cohomology 
problem has been solved using Morse theory ll| . These 
cases can also be solved using 'tic-tac-toe' with 5*1 and 

52 as in fig. [TJ To solve Hq^ , we start from the bottom 
chain. If a site on directly above this chain is occupied, 
there is an isolated site on the bottom chain, leading to a 
vanishing Hq^ . It follows that all sites directly above the 
bottom chain must be empty. Continuing this argument 
for subsequent chains, one finds that all sites on must 
be empty. Computing H12 is now a trivial step. The di- 
mension of Hq is related to the number of ground states 
on the chains that constitute 5*2. This depends on the ex- 
act boundary conditions. One finds that the dimension 
of Hq is cither one or zero, except for the cylindrical 
case periodic in the (to, — TO)-direction with to = 3p and 
n = 3g or n = 3(7 -|- 1. In that case the dimension of Hq 
is 2^, with K the nearest integer to 2n/3. 

For periodic boundary conditions in both diagonal di- 
rections, solving Hq becomes much more involved. This 
is due to the fact that in this case solving Hq^ does not 
imply that all sites on Si are empty. Instead, there are 
many allowed configurations on 5*1, and solving H12 be- 
comes highly non-trivial. We tackled this problem in two 
steps. First, we solved H12 for the case where S2 con- 
sists of 1 chain of arbitrary length, that is u = (to, —to) 
and V = (1,2). Then, we extended this to the case 
where 6*2 consists of an arbitrary number of chains, that 
is vi + V2 = 3p. This establishes ^ with A as given 
above. Our proof, which will be presented elsewhere, re- 
solves the ground state counting problem of this highly 
frustrated system on the 2D square lattice. Counting 
rhombus tilings, which is relatively easy 0, shows that 
indeed the number of ground states grows exponentially 
with the linear dimensions of the system. 

For V = (1,2), the tilings reduce to a linear sequence 
of tiles in arbitrary order. There are zero-energy ground 
states at all rational filling fractions in the range between 
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1/5 and 1/4. For example, for a system with 30 sites 
there are 5 Hnearly independent ground states with 6 par- 
ticles and 92 with 7 particles. The tilings not only count, 
but indeed seem to dominate the actual ground states. 
For specific lengths of the S2 chain, we have seen numeri- 
cally that the charge distribution of a ground state largely 
overlaps with that of the corresponding tiling. We can 
exploit this to gain physical insight. There are two uni- 
form phases: all squares at 1/5 filling and all diamonds 
at 1/4 filling. One diamond in a phase with all squares 
is then a zero-energy defect with fractional charge 1/5. 
From counting the number of tilings with one such defect 
it follows that a defect can have any momentum. This is 
reminiscent of a fiat band. Flat bands usually arise either 
from tuning the hopping terms of non-interacting parti- 
cles on an exotic lattice, or from tuning potential terms 
for strongly interacting particles with negligible kinetic 
energy [l^]. Here, however, the fiat band arises from 
tuning the potentials for particles with kinetic and po- 
tential energies of comparable size. Note that filling this 
flat band with defects is slightly subtle. Since defects 
cannot sit on top of each other and span over 4 sites, 
they obey a certain exclusion statistics. 
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FIG. 2: Edge modes. 



We can exploit the effective geometric description of 
the space of ground states even more by comparing pe- 
riodic and free boundary conditions. For the square 
lattice wrapped around the torus with u = (m, — m) 
and V = (n, n) the number of ground states grows as 
22(ra+m)/3 ^ |.]^g Cylinder, however, if one cuts the 
torus open along the w-direction only 2^"/^ ground states 
remain. Finally, if one also cuts the cylinder open along 
the w-direction one is left with a unique ground state on 
the plane. What happens to this vast number of ground 
states that disappear upon changing the boundary con- 
ditions? Consider the picture in fig. [21 If one identifies 
the dotted, zigzagged boundaries both in the u as well 
as in the w-direction, one finds that both the tiling with 
the drawn lines as well as the one with the dashed lines 
represent ground states. However, if one only identifies 
the boundaries in the M-direction, then the tiling with the 
drawn lines no longer represents a ground state. Instead, 
it has two defects at the edges, which can propagate along 



the edge. The only available scale for the energy of the 
edge mode is one over the length of the edge, which sug- 
gests that the edge modes are gapless. 

To further investigate criticality in this system, we 
have studied the model numerically with periodic bound- 
ary conditions such that u = {L,0) and v equals (0,2), 
(1, 2) or (3, 3). In all three cases, we have found com- 
pelling evidence for critical modes. We investigated how 
the spectrum changes upon twisting the boundary con- 
dition along the (L, 0)-direction from periodic to anti- 
periodic. This is a powerful way of distinguishing be- 
tween critical and gapped states: for a gapped state, the 
correlation length is finite and a change in the bound- 
ary conditions will have an exponentially small effect on 
the energy. In contrast, for a critical state the change in 
the energy will be substantial since the correlation length 
goes to infinity. More specifically, we will see that the en- 
ergy has a parabolic dependence on the boundary twist. 

In addition to this, one can extract more quantitative 
properties from twisting the boundary condition. For a 
critical supersymmetric system in ID with a Fermi sur- 
face, we expect that its continuum limit is described by 
an Af=2 superconformal field theory (SCFT). In such a 
theory, twisting the boundary condition corresponds to 
going from the Ramond to the Nevcu-Schwarz sector. 
The twist can be carried out continuously and leads to a 
spectral flow [l3|. If we define the twist parameter a to 
be integer in the Ramond sector and half-integer in the 
NS sector, the energy is a parabolic function of a, 



aQc 



(2) 



where c is the central charge. Qa depends linearly on a 
and is the sum of the left- and right-moving U(l) charges 
ll,[l3. The ir difference is conserved under the twist and 
is related to the fermion number. In the lattice model, 
we can go from periodic to anti-periodic boundary con- 
ditions continuously by replacing the term that hops a 
particle over the boundary cjyCi-l- h.c. by e^'^'^clyCi-t- 
h.c. The eigenvalues of the translation operator pa will 
now depend linearly on the twist parameter: 



^^), (3) 



so Pa = Po + oiF/ L where L is the length of the system 
and F is the total number of particles in the state \'4)). 

We computed the spectrum for various values of the 
twist parameter via exact diagonalization. We find that 
the majority of states have a parabolic dependence on the 
twist parameter by fitting a parabola to the energy levels 
as a function of the twist parameter, or equivalently, as 
a function of Pa ■ This clearly indicates that the system 
is critical. An example is shown in fig. [3l 

For a critical system, the energy of the SCFT in ^ is 
related to the numerically obtained value of the energy 
via -Bnum = '^t^EqftVf / L, where vp is the Fermi velocity 
and L the system size. So by comparing the parabolic 
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Momentum 

FIG. 3: In the plot we show nine energy spectra as a func- 
tion of momentum for increasing values of the twist parameter 
(0 < a < 1/2, with steps of 1/16) for a system with 33 sites 
(u = (11,0) and v = (3,3)) and 8 fermions. The lines are 
parabolic fits to the numerical data. The results from the 
spectral flow analysis for the red line through the origin ap- 
pear in red in table IT] 

We extracted values for E^/c and Qa/c via the above 
described method for three models with up to 36 sites 
(see table m. For the lowest energy levels, we typically 
find that (£'/c, Q/c) in the NS sector is either (-1/12, 0), 
(1/12, 1/3) or (1/4,2/3), all with an accuracy of within 
10%. These values occur in the Kac table for the fc-th 
minimal model of an M=2 SCFT with k even [li|. This 
is very compelling evidence that each of these systems is 
quantum critical. 

The fits become less reliable for levels with higher en- 
ergies, but also if there is an avoided level crossing as 
a function of the twist. This happens when the energy 
levels in the Ramond or NS sector are degenerate. For 
the chain, the avoided crossings vanish in the continuum 
limit 01 (see also [Hj), so one would expect that for the 
other models this is also merely a finite size feature. For 
the square ladder [u = (L,0) and v — (0,2)]. however, 
the results from exact diagonalization suggest that the 
avoided crossing will prevail for large system sizes. We 
are investigating this issue using density matrix renor- 
malization group methods [l^ . 

For the zigzag ladder [u = (L,0) and v — (1,2), see 
fig. 3], we have computed the entanglement entropy of 
the zero-energy ground states of the system up to 35 
sites. We find that for 1/4 and 1/5 filling, the entangle- 
ment entropy of the ground states saturates. This implies 
that the correlation length is finite, and thus the ground 
states at 1/4 and 1/5 filling are not critical. However, at 
intermediate fillings we find good correspondence to the 
spectral flow behavior of an M=2 SCFT. As anticipated 
in the introduction, this indicates that the system is crit- 
ical at 2/9 filling, where there is a phase transition due 



TABLE I: Results from spectral flow analysis for three types 
of tori (L, 0) X V. Here, A'^ denotes the number of sites and F 
the number of fermions. We show the results for the lowest 
energy level for each system. The values for E and Q are 
given in the NS sector {a = 1/2) and c is the central charge. 
For the level in red boldface the results are extracted from 
the red lino in fig. O 
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F jc 


Q/c 


18 


(3, 3) 
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-0.0851 


0.004 


36 


(3, 3) 
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-0.0841 


-0.002 


15 


(3, 3) 
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0.0898 


0.349 


21 


(3, 3) 
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0.0850 


0.337 


24 


(3, 3) 
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0.0850 


0.337 


30 


(3, 3) 


7 


0.0853 


0.338 


33 


(3,3) 
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0.0855 


0.338 


9 


(1.2) 


2 


-0.0858 


-0.005 


18 


(1,2) 
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r\ r\Ci Ac\ 

-0.0842 


-0.002 


27 


(1.2) 
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-0.0839 


-0.001 


17 


(1.2) 
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0.0844 


0.336 


26 


(1.2) 
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r\ r\Ci A r\ 

0.0840 


0.335 


35 


(1.2) 
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0.0839 


0.335 


14 


(1.2) 
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0.2666 


0.701 


23 




5 


0.2458 


0.657 


32 


(1.2) 
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0.2432 


0.652 


16 


(0.2) 
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-0.0897 


-0.014 


24 


(0.2) 
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-0.0889 


-0.012 


32 


(0,2) 
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-0.0885 


-0.011 


12 


(0,2) 
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0.0911 


0.350 


20 


(0,2) 
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0.0900 


0.348 


28 


(0, 2) 
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0.0894 


0.347 


14 


(0, 2) 
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0.0855 


0.338 


22 


(0,2) 
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0.0849 


0.337 


30 


(0,2) 


8 


0.0847 


0.336 



to competing orders of the diamonds and squares. 




FIG. 4: The figure shows the zigzag ladder [u = (i,0) and 
V = (1, 2)], where the zigzagged line is dashed to show that 
it can also be thought of as a chain with nearest and next- 
nearest neighbors excluded. 

The systems we have investigated numerically are 
rather confined in one direction. This suggests that we 
have essentially probed the edge modes and confirmed 
that they are gapless. Whether there are also gapless 
modes in the bulk remains unclear from this analysis. 
However, for a more generic case such as the triangu- 
lar lattice, where the ground state entropy is truly ex- 
tensive, one could speculate that a similar argument as 
above would imply that there are critical modes in the 
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bulk as well. 

It will be most interesting to investigate whether or 
not our main findings (quantum charge frustration and 
quantum criticality induced by strong repulsive interac- 
tions for itinerant lattice fermions) carry over to more 
generic models, and if they can be linked to some of the 
poorly understood physical features, transport proper- 
ties in particular, of strange metals and heavy fermions 
systems displaying non-Fermi liquid behavior. 
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